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ABSTRACT 

An update of a previous description of the BRUSLIB + NACRE package of nuclear data for astrophysics and of the web-based 
nuclear network generator NETGEN is presented. 

The new version of BRUSLIB contains the latest predictions of a wide variety of nuclear data based on the most recent version of the 
Brussels-Montreal Skyrme-Hartree-Fock-Bogoliubov model. The nuclear masses, radii, spin/parities, deformations, single-particle 
schemes, matter densities, nuclear level densities, EI strength functions, fission properties, and partition functions are provided for all 
nuclei lying between the proton and neutron drip lines over the 8 < Z < 110 range, whose evaluation is based on a unique microscopic 
model that ensures a good compromise between accuracy, reliability, and feasibility. In addition, these various ingredients are used to 
calculate about 100000 Hauser-Feshbach neutron-, proton-, a-, and y-induced reaction rates based on the reaction code TALYS. 
NACRE is superseded by the NACRE II compilation for 15 charged-particle transfer reactions and 19 charged-particle radiative cap- 
tures on stable targets with mass numbers A < 16. NACRE II features the inclusion of experimental data made available after the 
publication of NACRE in 1999 and up to 201 1. In addition, the extrapolation of the available data to the very low energies of astro- 
physical relevance is improved through the systematic use of phenomenological potential models. Uncertainties in the rates are also 
evaluated on this basis. 

Finally, the latest release vlO.O of the web-based tool NETGEN is presented. In addition to the data already used in the previous 
NETGEN package, it contains in a fully documented form the new BRUSLIB and NACRE II data, as well as new experiment-based 
radiative neutron capture cross sections. 

The full new versions of BRUSLIB, NACRE II, and NETGEN are available electronically from the nuclear database at 
pSttp ://w WW. astro. ulb. ac . b e/NuclearData| The nuclear material is presented in an extended tabular form complemented with a vari- 
ety of graphical interfaces. 
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■ 1. Introduction 

■ Through stellar evolution and stellar or non-stellar nucleosyn- 

■ thesis models, nuclear astrophysics is in high demand of a huge 
, variety of nuclear data. These include the static properties of nu- 
" clei, their spontaneous decay or fission characteristics, and their 
[ interactions with a variety of particles and nuclei. In particular, 

■ reactions of interest most often concern neutron-, proton-, or a- 
' particle induced reactions. Despite of significant efforts, experi- 
mental information only covers a minute fraction of the entire re- 
quired data. It is either insufficient, or even totally non-existent. 
Massive recourse to theoretical predictions is thus mandatory in 
many applications. 

The availability of compilations that provide an easy ac- 
cess to evaluated and well-documented nuclear data is an es- 
sential tool for astrophysics modelling, as is most dramatically 
illustrated by the impact on this type of research of the pio- 
neering compilation of charged-particle induced reactions by 
iFowler et al.l (Il967h . This has been the key motivation for the 
development of a new generation of compilations superseding 
[Fqwler et al. (1967) and its direct updates (Caughlan & Fowler 
(119881) . and references therein). The first such second genera- 
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tion compilation is referred to as NACRE (Nuclear Astrophysics 
Compilation of REactions) (lAngulo et al. It includes 

several new important features with respect to its predeces- 
sors. As an extension of NACRE, a library of nuclear data re- 
ferred to as BR USLIB (BRUSsels LIBrary) has been constructed 
dAikawa et al.L .2005) . On top of the NACRE data, it includes a 
wide variety of nuclear information of interest, and takes full 
advantage of the internet facility, which allows one to make 
a huge volume of data available door-to-door in an accurate 
and usable format. In addition, a nuclear NETw ork GENerator 
NETG EN complements the BRUSLIB package dAikawa et all 
l2005h . Recently, the Karlsruhe Astr ophysical Database of 
Nucleosynthesis in Stars (KADoNiS) (iDillmann et al.' '2006a!) 
and the JINA REACLIB database (Cyburt et al., 2010) have 
been constructed and are updated regularly. 

The aim of this paper is to present an update of BRUSLIB 
and NETGEN. BRUSLIB now contains new predictions based 
on microscopic models. As presented in Sect. 2, these include 
nuclear masses, deformations, single-particle schemes, density 
distributions, nuclear level densities, partition functions. El 
strength functions, and fission properties. These data are used 
to calculate a large body of thermonuclear reaction rates based 
on the Hauser-Feshbach (HE) model. BRUSLIB also contains 
an update of NACRE, referred to as NACRE II, that is now com- 
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Fig. 1. (Colour online) HFB-21 nuclear mass excesses. 
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pleted dXuet all l2009l 1201 2|) . As discussed in Sect. 3, NACRE 
II revises 34 NACRE rates based on experimental data pub- 
lished after the NACRE release, and on a new theoretical han- 
dling of these data. Section 4 is devoted to the description of the 
NETGEN update. A summary and some prospects are presented 
in Sect. 5. 



2. BRUSLIB 

2.1. Nuclear ground-state properties 

Microscopic models based on realistic nucleon interactions have 
been used to estimate various nuclear ground-state properties. 
One of the most promising approaches of this sort is the non- 
relativistic Hartree-Fock-Bogoliubov (HFB) method, which is 
based on an effective nucleon-nucleon interaction of the Skyrme 
type. HFB calculations in which a Skyrme force is fitted to es- 
sentially all the mass data dAudU I2OIQ) can now successfully 
compete with the most accurate droplet-like formula. 

The Skyrme HFB mass model not only seeks optimized fits 
to the mass data, but it also allows the construction of a universal 
effective interaction capable of reproducing most of the observ- 
ables of relevance in nuclear applications. This essential prop- 
erty is obtained by imposing extra physical constraints on the 
HFB mass model. In its last version, additional terms are inserted 
in the effective nucleon-nucleon interaction to better account for 
different physical effects, among which the pairing interaction, 
or the properties of neutron matter for the description of the in- 
ner crust and the core of neutron stars. Details can be found in 



Tondeur et al. (2000), Samvn et al. (2001), Gorielv et al. (2002), 
Gorielv et al._ (2003) Samvn et al (.20041 Gorielv et al, (2005.) 



Gorielv et al 



Gorielv et al 



(feooi.lGodelvetalJ(^^ 



2009ah . and lGorielv et al 



JChamel et al.l(l2008h . 
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So far, 21 HFB mass tables have been constructed with 
different par a meter izations of the Skyrme effective force. 
iGorielv et al.l (l2010h presented the most recent HFB-21 mass 
values and related ground-state properties. Tables of 8508 
masses containing all nuclides with 8 < Z < 110 lying between 
the proton and the neutron drip lines derived from all 21 HFB 
models can be accessed from BRUSLIB. For each nuclide, the 
calculated mass is complemented with its nuclear deformation 



parameters y32 and y64, rms charge radius R^h, deformation energy 
Edef, neutron S„ and proton S p separation energies, beta-decay 
energy Qhei, the difference Mgrr between the experimental and 
calculated mass excesses as well as the experimental Jt,xp and 
calculated J,/, ground-state spins, and the experimental P^^p and 
calculated f ground-state parities. Figure [T] shows the latest 
HFB-21 mass excesses in the nuclear chart. 

By solving the HFB equation, the single-particle spectra, nu- 
clear potential, and charge- and matter density distributions can 
also be obtained for each nucleus. BRUSLIB provides these re- 
sults for all 8508 nuclides. These basic nuclear properties pro- 
vide the essential input required for additional predictions. As 
examples, the nuclear density distribution is used for micro- 
scopic or semi-microscopic optical model potential calculations 
by folding the target radial matter density with an optical poten- 
tial in nuclear matter, and the single-particle levels are employed 
for evaluating the nuclear level density (see Sect. 2.2). 



2.2. Nuclear level densities and partition functions 

The nuclear level density (NLD) is an essential ingredient for 
cross section and reaction rate calculations within the statistical 
HE model. A microscopic combinato rial approach for estimating 
the NLDs has been first described bv lHilaire et"aD (|l998, 200 1|. 
It has been improved subsequently to include both collective ef- 
fects and pairing correlations (Hilaire & Gorielv. 2006). 

The model uses the constrained axially symmetric 
HFB method based on the BSkl4 effective Skyrme force 
(iGorielv et al.L l2007h to construct incoherent particle-hole state 
densities as functions of the excitation energy, spin projec- 
tion, and parity. In its latest developments, the combinatorial 
method improves the description of the collective vibration 
levels by taking the phonon excitation explicitly into account 
(Goriely et al., 2008b). 

BRUSLI B provides th e spin- and parity-dependent NLDs 
from Gorielv et all (l2008bl) for 8508 nuclides with 8 < Z < 1 10 
lying between the proton and the neutron drip lines, for excita- 
tion energies U up to 200 MeV, and for spins up to J = 49 for 
even-A nuclei or 99/2 for odd-A nuclei. The nuclear temperature, 
cumulative and total numbers of levels are also included. No 
simple analytical fits to the tabulated NLDs are given to avoid 
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Fig. 2. (Colour online) Spin- and parity-dependent NLDs for ^''Fe, ''''Ru, ^'^Sm and ^'^Bi predicted bv lGorielv etal] (l2008bl) . 



losing the specific microscopic characteristics of the model. A 
re-normalization procedure of the NLD on experimental data is 
required in many instances, and in particular for nuclear data 
evaluation or for an accurate and reliable estimate of reaction 
cross sections. More specifically, the re-normalized level densi- 
ties are derived through the formula 



p(U, J") ^p(U-6, J") exp(a Vt/ - 6), 



(1) 



where the energy shift 6 of the excitation energy U is essen- 
tially extracted from the analysis of the cumulative number of 
experimental levels, and a from the experimental s-wave neutron 
spacings. Equation ([TJ has been used to fit the 289 nuclides for 
which this experimental information exists. The corresponding 6 
and a values along with the experimental data are also given in 
BRUSLIB. Figure |2]illusti-ates the calculated ^^Fe, ^^Ru, "^'^Sm, 
and ^''^Bi NLDs for different spins and parities. For these same 
nuclides, Fig.[3]shows that the predicted cumulative numbers of 
levels reproduce the measured values quite well. 

The astrophysical reaction rates require the knowledge of the 
temperature-dependent partition function for each nuclide (Z, A). 
Its value normalized to the ground-state value is given by 



GiT) = 2 



2Ji+l ^ U: 

exp( ) 

2Jq + 1 ^ kT 



/ 



^ 2J+ I U 



2Jo + V 



(2) 



where Jq is the ground-state spin, and f/, and 7, are the excitation 
energy and spin for the \-th discrete excited state. The summation 
extends over the experimentally known energy levels up to an 
excitation energy t/,/,, above which the knowledge of the energy 
spectrum is considered to be incomplete. The integration covers 
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Fig. 3. (Colour online) Comparison of the cumulative num- 
bers of ener gy lev els derived from experiment (solid lines) 



(Capote et al.Ll2009h and those predicted (dashed line) by the mi- 
croscopic NLD model of lOorielv et al.. aOOSb) for ^^Fe, '''^Ru, 
i44Sm, and ^o^Bi. 



the energy range above f/,/„ and involves levels with all possi- 
ble spins J of parities ji - + and -. In practice the BRUSLIB 
partition functions are calculat ed using the experirn ental excited 
spectrum whenever available ICapoteeTaD, 120091) and, other- 
wise, the NLDs predicted by Goriely et al. (2008b). BRUSLIB 
provides the partition functions for the same set of 8508 nu- 
clides, and for temperatures in the 0.001 < < 10 range (where 
Tt) denotes the temperature in lO'^K). For illustration. Fig. |4] 
shows the partition functions for ^*Fe, ^''Ru, '"^"^Sm, and ^"^Bi. 
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Fig. 4. (Colour online) Partition functions for ^''Fe, ^^Ru, '^''Sm, 
and ^°^Bi as a function of the temperature. 



2.3. E1 strength function 

Gamma-ray strength functions are involved in the HF calcula- 
tion of capture cross sections, y-ray production spectra, isomeric 
state populations, and competition between y-ray and particle 
emission. In this context, relevant multi -polarities are El, Ml, 
and E2 strengths, among which El dominates in general. 

Large-scal e derivations of the El strength function 

dGorielv et all 12004) have been conducted with the use of the 
HFB pl us quasi-par t icle ra ndom phase approximation (QRPA) 
models (iKhan et a n 1200 Ih based on a realistic Skyrme inter- 
action. The HFB model allows one to treat pairing effects on 
the ground state in a self-consistent way, while QRPA views 
the collective nuclear excitation as a collective superposition of 
two quasi-particle states built on top of the HFB ground state. 
This collective aspect of the excitation makes QRPA an accurate 
tool for investigati ng the El -stre ngth function in both closed and 
open shell nuclei dGorielv et all 2004). BRUSLIB contains the 
El strength functions obtained from the HFB plus QRPA mod- 
els based on the BSkl4 Skyrme force for all 8508 nuclides with 
8 < Z < 1 10 lying between the two drip lines. Data are presented 
on an energy grid of 0. 1 MeV between and 30 Me V. 

Figure |5] illustrates the El-strength functions for some Sn 
isotopes. In the neutron-deficient region as well as along the 
jS-valley of stability, the strengths are very similar to the em- 
pirical Lorentzian-like approximation. In the neutron-rich re- 
gion, the HFB-i-QRPA predictions start deviating from a simple 
Lorentzian shape with some extra strength located at an energy 
lower than the giant dipole resonance energy. The more exotic 
the nucleus, the stronger this low-energy component referred to 
as the "El pygmy resona nce". This prediction agrees with exist- 
ing experimental results dGorielv et al.Ll20()4h . 

2.4. Fission properties 

Fission is of special importance in the modelling of the r-process 
of nucleosynthesis, which is invoked to explain the origin of 
approximately half of the stable nuclides above iron observed 
in nature, as well as the production of the long-lived actinides 
^^^Th, ^^^U, and ^^**U, which are often used as cosmochronome- 
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Fig. 5. (Colour online) El-strength functions of some Sn iso- 
topes. 

ters to estimate an upper limit of the age of the Galaxy. For 
this purpose, the probabilities of spontaneous, j6-delayed and 
neutron-induced fission for the nuclides with 80 < Z < 110 
mainly located in the neutron-rich region of the nuclear chart are 
needed. 

All key ingredients entering the fission description, such as 
the static fission path (including the height and width of the pri- 
mary, secondary and possible tertiary barriers) and the NLD, are 
derived within the same microscopic model. The HFB model 
with the BSkl4 Skyrme force is used to estimate the static fis- 
sion path (allowing for axially symmetric deformations as well 
as left-right asymmetries) for all nuclides with 90 < Z < 110 
lying between the vall ey of j6-stability and the neutron drip lines 
dGorielv et al.L l2009bh . These predict ions have been show n to 
agree with existing experimental data dGorielv et al.Ll20()9 b'). 

Figure |6] displays the fission paths for the U and Cm iso- 
topes close to the valley of /? stability. Each fission path corre- 
sponds to the most gently climbing or steepest descending path 
found and projected along one deformation parameter, namely 
the quadrupole deformation /32. Although a tiny third barrier 
clearly appears at strong deformations for some U isotopes, the 
fission path for these nuclei appears to be well-represented by 
a traditional double-humped barrier, at least in a local region 
located close to the saddle-point deformations. The situation 
can be quite diff'eren t farthe r away from stability dGorielv et al.L 
l2008d: ICapote et all l2009h . In fact, the fission path for exotic 
neutron-rich nuclei cannot be simply approximated by a double- 
humped barrier with parabolic shapes. 

The HFB plus combinatorial method developed to estimate 
the NLD at ground-state deformation (Sect. 2.2) is also used to 
calculate the NLD at the saddle points along the fission path, 
making a coherent use of the corresponding HFB predictions for 
the single-particle level scheme and pairing strength at the cor- 
responding deformation (Gorielv et al., 2009b). Because of the 
lack of experimental information, the same prescription is used 
for the saddle points as for the ground state, i.e., a total of three 
phonons are coupled to the excitation configurations of a maxi- 
mum of four particle-holes. Quadrupole, octupole, and hexade- 
capole phonons are included, their energies being assumed iden- 
tical to those of the ground state. BRUSLIB provides the corre- 
sponding NLDs for 90 < Z < 110 nuclides at each of their two 
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Fig. 6. (Colour online) Predicted HFB fission paths for the Cm Fig. 7. (Colour online) Total NLDs for ^^^Cm and ^^^V at diff'er- 
and U isotopes close to the valley of/3 stability versus the defor- ent locations along their fission path, 
mation parameter /32- 



(or three) highest saddle point barriers and one (or two) shape 
isomers, and at energies up to f/ = 200 MeV and spins up to 
7 = 49 (99/2). As an example, the ^*^Cm and ^^^U total NLDs 
are shown in Fig.Qat different locations along the fission paths. 



3. HF reaction rates 

Thousands of reactions of astrophysics relevance involve excited 
states of the target, and more or less exotic nuclides. No wonder 
then that experimental data are largely missing, and will remain 
so for a long time to come. Theory has therefore to provide the 
necessary data. 

Th e reaction rates supplied by BRUSLIB dGorielv et al.L 
I2008ah are obtained with the code TALYS, which is ded- 
icated to nuclear reaction simulations through the use of 
many state-of-the-art nuclear models ( Koning et all l2002t 
iKoning & DuiivestiinL l2004t iKoning et al.L 12004'). It is able to 
cover all main mechanisms at play in light particle-induced nu- 
clear reactions, and provides a complete description of all re- 
action channels. In particular, TALYS takes into account most 
types of direct, pre-equilibrium, and compound mechanisms to 
estimate the total reaction probability, as well as the competi- 
tion between the various open channels. The code is optimized 
for incident projectile energies ranging from 1 keV to 200 MeV 
on target nuclei with mass numbers between 12 and 339. It in- 



cludes photon, neutron, proton, deuteron, triton, -^He, and a- 
particles as both projectiles and ejectiles, and single-particle as 
well as multi-particle emissions and fission. All experimental in- 
formation on nuclear masses, deformations, and low-ly i ng sta tes 
spectra from, e.g., the RIPL-3 database (Capote et al., 2009) is 
adopted as input. If the necessary laboratory information is miss- 
ing, the BRUSLIB data described in Sects. 2. 1-2.4 are used. This 
concerns the nuclear masses, spins, parities, NLDs at ground 
state or fission saddle points, and partition functions. In addi- 
tion, the El strengths (Sect. 2.3) were used to calculate the pho- 
ton transmission coefficients Ty, and the fission path s of Sect . 
2.4 were used in conjunctio n with the WKB method (ISin et al.L 
l2006HGorielv et al.Ll2009bl) to calculate the fission transmission 
coefficients Tf. As for the particle emission channels, the trans- 
mission co efficients were eva luated on the basis of the optical 
potential of lKoning & Delarocha (.2002.) fo r the nucleon-nucleus 
interaction and of Demetriou et al. r il200l for the a-nucleus in- 
teraction. 

BRUSLIB supplies about 100000 TALYS rates for (n,y), 
(n,p), (n,^), (p,r), (p,n), (p,a), ia,y), (a,n), ia,p), (y,n), (y,p), 
and {y,a) reactions on some 8500 targets with 8 < Z < 110 
lying between the proton and the neutron drip lines, and for 
temperatures in the 0.001 < Tg < 10 range. For nuclides with 
80 < Z < 110 that lie between the proton and the neutron drip 
lines, each of the neutron-, proton-, a-, and y-induced fission 
rates, i.e. the channels (n,f), (p,f), (a,f), and (y,f), can be ac- 



6 



Y. Xu, S. Goriely, A. Jorissen, et al.: BRUSLIB, NACRE 11 and NETGEN 




cessed from BRUSLIB. For Z > 90, they are calculated on the 
basis of the HF model making use of the fission barriers and the 
ground-state and saddle point NLDs calculated as described in 
Sect. 2.4. As an example, Fig.[8]shows the various reaction rates 
on the targets '^°Sn and ^^^U. 

The TALYS estimates for n-, p-, deuton-, triton-, ''He-, and 
a-induced reaction rates on targets from Li to Na that c ould 
be involved in Big-Bang nucleosynthesis dCoc et all 1201 Ih are 
also included in BRUSLIB. These predictions can be viewed at 
best as rough first approximations when no information from 
other sources is available. As summarized bv lCoc et al.l (1201 ll) . 
TALYS is expected to provide these reaction rate predictions 
within a 3 orders of magnitude accuracy in the temperature range 
of Big-Bang nucleosynthesis, though in many cases the agree- 
ment can be regarded as satisfactory. 

More generally, the uncertainties in the calculated cross sec- 
tions arise from two main origins. The first one is related to 
the reaction mechanism itself. All codes previously developed 
for astrophysical applications only describe one reaction mecha- 
nism, either of direct or of compound nucleus type, and include a 
number of approximations. Even if many of these can be justified 
by the very low relative energies (often a few keV to 1 MeV) of 
astrophysical relevance, it is of interest to test their validity on a 
quantitative basi s. TALYS is suited for this, in particular through 
the inclusion of (iGorielv et al.L l2008ah 

- the pre-equilibrium reaction mechanism; 

- the detailed description of the decay scheme, including the 
description of y-delayed particle emission and the possible 
particle emission from all residual nuclei; 

- the multi-particle emission; 

- the choice of a variety of width-fluctuation correction fac- 
tors; 

- the coupled channel description for deformed nuclei; 



- the fission channel for the compound as well as for the resid- 
ual nuclei. 

Other uncertainties come from the nuclear ingredients that en- 
ter the reaction model. When dealing with nuclear astrophysics 
applications, these ingredients should be derived from global, 
universal, and microscopic models. The large number of nu- 
clei involved in modelling of some nucleosynthesis mecha- 
nisms indeed implies that global models should be used. On 
the other hand, a universal description of all nuclear proper- 
ties within a unique framework for all nuclides involved in 
a nuclear network ensures the essential coherence of the pre- 
dictions for all unknown data. Finally, a microscopic descrip- 
tion provided by a physically sound theory based on first prin- 
ciples likely renders extrapolations away from experimentally 
known energy or mass regions more reliable than predictions 
derived from more or less parametric approaches of various 
types and levels of sophistication. The microscopic models used 
in TALYS to calculate the ingredients of the BRUSLIB reac- 
tion rates meet these criteria through a satisfactory compro- 
mise between accuracy and reliability. This would make the 
BRUSLIB reaction rate library unique away from stability, while 
it is just of comparable reliability close to stability, as illus- 
ti-ate d by Fig. H The f uU BRUSLIB database can be found at 
http: //www. astro, ulb. ac. b e/bruslib 

4. Experiment-based reaction rates: the NACRE II 
compilation 

Since the NACRE compilation dAngulo et all Il999h . many 
charged-particle-induced cross sections of astrophysical interest 
have been measured or re-measured. In addition, it has been con- 
sidered desirable to invest more effort than into NACRE on the 
modelling of the extrapolations of the experimental data to the 
very low energies of astrophysical relevance. These extrapola- 
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Fig. 9. (Colour online) Comparison between the T<) - 0.3 
BRUSLIB radiative neutron capture rates and t he ex perimental 
rates provided by the compilation of iBao et"al] (1200 0). comple- 
mented with data published in refereed journals between 2000 
to 201 1 (references mentioned in Sect. 5). 

Table 1. Reactions included in NACRE II. 



No. 


reaction 


No. 


reaction 


1 


'H(p,y)'He 


18 


''Be(p,7)'"B 


2 


^U(d,yflle 


19 


''Be(p,d)2''He 


3 


2H(d,n)3He 


20 


''Be(p,a)«Li 


4 


2H(d,p)3H 


21 


''Be(Q',n)'2c 


5 


^H(a,yfU 


22 


"'B(p,7)"C 


6 


3H(d,n)''He 


23 


"'B(p,Q')'Be 


7 


3H(a,7)'Li 


24 


"B(p,7)'^C 


8 


3He(d,p)'*He 


25 


"B(p,Q-)2-*He 


9 


3He(3He,2p)'*He 


26 


"B(c!',n)"'N 


10 


3He(Q',7)^Be 


27 


'-C(p,7)"N 


11 


*Li(p,7)'Be 


28 


'2C(Qr,7)"*0 


12 


'*Li(p,a)3He 


29 


'3C(p,7)'-*N 


13 


^Li(p,7)'*Be 


30 


'3C(Q',n)'«0 


14 


'Li(p,a)'*He 


31 


"N(p,7)''*0 


15 


^Li(a,7)"B 


32 


'*N(p,7)'50 


16 


■'Be(p,7)«B 


33 


"N(p,7)'*0 


17 


'Be(a',7)"C 


34 


'5N(p,q')'2c 



tions are essential because the relevant energy domain cannot 
be reached in the laboratory, except in very few cases. This sit- 
uation has motivated the launch of the NACRE II compilation, 
which aims at re v ising some of the NA CRE rates (Aikawa et aL, 
2006; Katsumal l2006t lArnould & KatsumaL 120071: iXuetal. . 
2009.2012) . 

In its current version, NACRE II surveys 34 two-body 
exothermic reactions (15 particle-transfer and 19 radiative cap- 
ture reactions) on targets with A < 16 listed in Table [1] 
Phenomenological potential models are adopted to describe and 
extrapolate their resonant and non-resonant cross sections at 
very low energies. The HE approximation is not suited for the 
cases under consideration, except possibly at the highest consid- 
ered energies (or temperatures). 

The methodology adopted for the construction of NACRE II 
can be summarized as follows: 



Reaction mechanisms Two main mechanisms are considered for 
low-energy nuclear reactions: the compound nucleus and the di- 
rect reaction processes. The direct reaction process is important, 
and often dominant, in charged-particle-induced reactions at the 
very low energies of astrophysical interest. Given the difficulty to 
tunnel through the Coulomb barrier, the reaction may occur be- 
fore the projectile can penetrate deep inside the target nucleus. 
The formation of a compound state is accordingly suppressed, 
which is reinforced especially in light nuclei by the paucity of 
quasi-bound levels. 

Adopted models The potential model, including the El, E2, and 
Ml transitions, is used to evaluate the 19 capture reactions, while 
the DWBA with zero-range interaction approximation is em- 
ployed to study the 15 transfer reactions. The adopted nuclear 
potential is a real Woods-Saxon potential with a surface ab- 
sorption imaginary part. At low energies, the significant direct 
and various discrete resonant contributions to the cross sections 
are both taken into account in these two models. However, for 
a few cases involving extremely narrow resonances, the Breit- 
Wigner formula is adopted for evaluating the corresponding re- 
action rates around the resonance energy. 

Selection of experimental data The primary ensemble of exper- 
imental low-energy cross section data of NACRE II comprises 
those included in NACRE and supplementary ones that have 
been published by the end of 201 1 preferentially in refereed jour- 
nals. Generally speaking, we take the selected experimental data 
on cross sections and associated errors at face value since we are 
not confident enough to do otherwise, given the quite limited in- 
formation that is normally available to us. We also do not refer 
in NACRE II to any differential quantities, even when they have 
been measured. 

Fitting procedure For transfer reactions, a DWBA model is 
adopted with a cut-off energy for the free parameter fit on the 
high-energy side normally set at Ec.,,,. - 1 MeV. However, when 
resonances are found experimentally below 1 MeV, only the one 
with the lowest excitation energy is embraced into the fit. This 
procedure is justified because the S-factor below that energy 
range normally behaves like the tail of that resonance. Just for 
the practical purpose of reproducing the measured cross section 
data, one may drastically reduce the number of potential parame- 
ters of the entrance and exit channels and of the form factor with- 
out causing much damage. First of all, for the entrance channel, 
we generally retain only a shallow, imaginary part of the nu- 
clear potential, which takes into account the weak absorption to 
the exit channel by particle transfer We then parameterize rather 
than fit the radii and difFuseness of the global potentials. The po- 
tential depths are left as the adjustable parameters to be fitted. 
We also treat the spectroscopic factor as an adjustable parameter 
for the absolute value of the cross section. In general, the best 
values of the adjustable parameters have been derived by apply- 
ing the standard fit with the occasional help of a fit-by-eye. 

A potential model is used for the calculation of radiative cap- 
ture cross sections. These are rarely measured below Ec.,,,. =0.1 
MeV, where non-resonant contributions may well be dominant. 
Therefore we adopt a strategy of parameterization that is quite 
different from the one selected for transfer reactions. In partic- 
ular, we try to fit the resonances, if any, with the real potential 
form to deduce the non-resonant contributions simultaneously. 
To reduce the number of parameters, we take the same set of 
radii and difFuseness for both the initial and final states. Hence, 
for a given resonant or non-resonant contribution, the potential 
depth, radii, diffuseness, and the re-normalization constant in- 
cluding the spectroscopic factor are left to be fitted. This is quite 
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in contrast to the strategy used for transfer reactions, and fit-by- 
eye is used more often than in the case of transfer reactions. 
Reaction rate evaluation For each reaction, three different sets 
of model parameters are determined, duly taking into account re- 
ported experimental uncertainties. From these sets, the adopted, 
low, and high limits of the extrapolated S-factors are calculated 
by the reaction model. Reaction rates Na < crv > are evaluated 
in the 0.001 < < 10 temperature range. All details can be 
found at http://www.astro.ulb.ac.be/nacreii 

Many mor e details on NAC RE II than those provided above 
can be found in fXu et al.l ( l2012h . In particular, specific comments 
are given for each of the reactions of Table [1] along with the 
values of all model parameters used to calculate the tabulated 
adopted, lower and upper limits of the rates. The NACRE II rates 
are also compared to the NACRE rates. 

5. NETGEN 

5.1. Content 

The nuclear NETwork GENerator NETGEN vlO.O is an inter- 
active, web-based tool to help astrophysicists in building up 
a nuclear reaction network as defined by each user. It gener- 
ates tables of the necessary nuclear reaction rates on a tem- 
perature grid specified by the user, and provides the references 
of the sources of these rates. All reaction rates, theoretical or 
experiment-based, include the contribution of thermally popu- 
lated excited target states through the calculation of the stellar 
enhancement factor. 

The following nuclear data are available in NETGEN for 
about 8600 nuclides located between the neutron and proton drip 
lines in the 1 < Z < 110 range: 

- for compiled experiment-based nuclear reaction rates: 

- the NACRE II dat a (see Sect. 4) fori < A < 16 targets; 

- the NACRE data (l Angulo et al. l [T999l) . when not super- 
seded by NACRE II, or by other more recent compila- 
tions than NACRE listed below; 

- the compilation of the r ates of Big-Bang reactions b ased 
on the R-matrix model dPescouvemont et"ani2004l) : 

- the 2010 evaluation of Monte Carlo-based rates of reso- 
nant proto n capture re actions on targets in the 20 < A < 
40 range (Iliadi s etali l2010); 

- the proton-induced reaction rates on stable and unstable 
target nuclei in the 20 < A < 40 region (llhadis et al.L 
I2001h : 

- when not available in the above- mentioned compila- 
tions, the reaction rates proposed bv lCaughlan & Fowled 
(11988,) : 

- the compilation of iBao et"al ] (l2000h for radiative neutron 
capture rates; 

- various radiative neutron capture rates made available 
between 2000 and 2011 that superse de the c o mpila- 
tion of iBao et al.l (l2000h (TAbbondan no et all l2004| 
Aerts et al.L l2006t iBest et al.L l200ir :Blackmon et al., 
20021; 'Borella et al.', '2007; 'Dillmann et al., 2006b, 
2009, 2010; Domingo-Pardo et al., 2006a b, 2007a b; 
Ischetal., 2008; Guberetal., 2002, 2010; Heil et al., 
2005, 2008a b e; Ka toh et al..,20q3; Koehler et aL.1998) 



200^Reifarth et al., 2002^ 120081: ISedvshev et al.L 1 1999 : 



zvuj, zuuoa u c, rs^a ioii ei ai.».iz,uuj. rs^oeinei ei ai.^ iiyyoL 
2(300, 2001; Ledere ret^ 1201 U iMaiTone etail 2006[ 
Marganiec et al., 2009a b, 2010; Meierhofer et al.,7oH 



Shcherbakovetal., 2005"; 'Tagliente et al.', '2008a'b 
2010 2011; Terlizzietal., 2007; Uberseder et al] 
2009t (Voc kenhuberetal.. 2006; Voss et al., 1999 
Winckleii 12006: .Wisshak et al.u200ai200k, 2002. .2004 



2006allbllch 

additional specific experimental reaction rates published 
in refereed journals that were not included in any of the 
above compilations. 



for theory based nuclear reaction rates: the rates of the 
neutron-, proton-, a-, and y-induced reaction rates, includ- 
ing fission, provided by BRUSLIB (Sect. 2). 

for compiled experimental and theoretical y6-decay and elec- 
tron capture rates: 

- the experimental ;8- decay rates of Katakura J 19961) ; 

- the ,6-decay ra tes of lTakahashi & Yokoil ( 1987 ) (see also 
iGorielvl [19991) : 

- specific electron capture rates from various published 
sources; 

- ; 6-decav rates ca l culate d with the revised Gross theory 
dTachibana et all 1 19901) and Q- values from the HFB-14 
mass model; 

- jS-decay r ates calculated with the QRPA plus FRDM 
models bvl MoUer etal] (Il996h ; 

- the density- and temp erature-dep e ndent /3- 
decay rates calculated by Oda et al.l (1 19941) or 
iLanganke & Martinez-Pinedol(l2000h . 



Mphr et al., 1999; Mosconi et al., 2010; Nakamura et al., 
20031; iNo guereetal., 2006; O-Brien et al., 2003; 
Patronis et al.>. i2004t iRapp et al.t. .2002; ,Ratzel et al.L 



5.2. Instructions for use of NETGEN 

For each nuclear reaction, yS-decay and electron-capture rate on 
a nuclide involved in a network defined by the user, NETGEN 
adopts by default the data source that is considered to be the most 
reliable. In order of preference, it selects the latest available com- 
pilation based on experimental data, followed by the BRUSLIB 
rates. The user may nevertheless make another choice for se- 
lected cases by specifying a bibliographic index for each reac- 
tion. The possible choices for this bibliographic index are given 
in a help panel on the website and in the output log file. NaN 
(Not a Number) flags are present in the output if the tempera- 
ture grid requested by the user extends beyond the limit of valid- 
ity of the experimental range. To remove these NaN, NETGEN 
provides the option to use theory-based reaction rates properly 
scaled (to ensure continuity). Finally, a FORTRAN computer 
program handling the rates is also provided. More details can 
be found at http://www.astro.ulb.ac.be/Netgen. 

The general steps for constructing a network and the associ- 
ated reaction rate library are as follows: 

Step 1: Network preparation A network may be generated by 
four non-mutually exclusive methods: 

- reactions typed in one by one. When choosing a reaction, the 
bibliographic index should also be given; 

- automatically generated network in the (A^, Z) plane. The nu- 
clides involved can be selected by charge (Z) and mass (A) 
number, boundaries (drip lines, stability line), and among 
various possible sets (proton-, neutron-, (T-captures, photo- 
dissociations, etc.). In this option, the default reaction rates 
{bibliographic index =1) are selected. This means that if sev- 
eral rates are available for a given reaction, the experimental 
or experiment-based rate is chosen rather than the theoretical 



Y. Xu, S. Goriely, A. Jorissen, et al.: BRUSLIB, NACRE 11 and NETGEN 



9 



one. If several experimental rates are provided, the most re- 
cently published one is selected. For ,8-decays, the rates are 
taken from Takahashi & Yokoi (1987). If a rate is not avail- 
able from this source, other experimental and most recent 
rates are chosen; 

- upload the reaction file according to the NETGEN prescrip- 
tion; 

- select the reactions from the specific pre-defined database. 

Step 2: Network construction After step 1, the user has the pos- 
sibility to list the current network, modify it, sort it, or clear it. 
Plotting the network in the (A^, Z) plane and individual reaction 
rates is offered as additional options. 

Step 3: Compute the reaction rates The network reaction rates 
are computed on the user-defined temperature grid. A file record- 
ing the rates is generated, as well as a log file that provides the 
references from which the rates have been calculated. Plots of 
the rates versus temperature can be obtained. 

6. Summary and outlook 

The last update of the nuclear databases BRUSLIB 
and NACRE II, and of the nuclear network generator 
NETGEN was presented. This package is of unprecedented 
scope and quality. It is accessible through the website 
\http ://www. astro, ulb. ac. be/N uclea rData\ 

BRUSLIB contains an immense variety of static properties 
for nuclei that lie between the proton and neutron drip lines over 
the range 8 < Z < 110, as well as the rates of proton-, neutron- 
, and a-particle-induced reactions on these nulides. Calculated 
rates of interest for Big-Bang nucleosynthesis are also given, 
but have to be considered as rough guesses only. BRUSLIB has 
the unique feature that this vast body of data are predicted by 
a global and universal microscopic model, and by the modern 
reaction code TALYS. All BRUSLIB data can be accessed via 
user-interactive figures and tables. 

NACRE II supersedes NACRE for 15 transfer and 19 radia- 
tive capture reactions on nuclides with A < 16. It includes ex- 
perimental data published after the release of NACRE and up to 
2011, and makes use of a phenomenological model to extrap- 
olate the experimental data to the very low energies of astro- 
physical relevance, and to evaluate the uncertainties on the rates. 
Detailed electronic material is supplied, in particular, tabulated 
rates. 

The web-based nuclear network generator NETGEN gives 
the possibility to set up a network, and to obtain the necessary 
network (nuclear and weak) reaction rates in a tabular form on a 
grid selected by the user, and from sources that are selected by 
default, or by the user 

The databases described in this work of course need contin- 
uous updating and improvement. Just as an illustration, TALYS 
offers the possibility of using alternative local or global mod- 
els for predicting various nuclear ingredients. Additional ef- 
fort should be put in the study of not only parameter uncer- 
tainties, but also of model uncertainties (e.g. the a-particle 
optical potential) that affect the astrophysica l reaction rates 
dArnould & Goriel v','2006V Arnould et aUl2007l) . In addition, di- 
rect captures not treated in TALYS should be scrutinized on a 
large-scale basis with the help of a global model using of the 
same ingredients as those used in HE calculations. The release 
of an extended spallation reaction library is also envisioned. 
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